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The anisotropies in the pressure obtained from the energy-momentum tensor are studied for 
magnetized quark matter within the su(3) Nambu-Jona-Lasinio model for both /3-equilibrium matter 
and quark matter with equal quark chemical potentials. The effect of the magnetic field on the 
particle polarization, magnetization and quark matter constituents is discussed. It is shown that 
the onset of the s-quark after chiral symmetry restoration of the u and d-quarks gives rise to a 
special effect on the magnetization in the corresponding density range: a quite small magnetization 
just before the s onset is followed by a strong increase of this quantity as soon as the s quark sets in. 

It is also demonstrated that for B < 10 18 G within the two scenarios discussed, always considering 
a constant magnetic field, the two components of pressure are practically coincident. 
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The structure of the QCD phase diagram is of utmost importance in understanding many physical aspects of 
nature, ranging from the early universe to possible nuclear liquid-gas and hadronic quark matter phase transitions to 
the physics of compact objects [ljj. Early analyzes performed within the Nambu -Jona-Lasinio model (NJL) framework 
indicate that when strongly interacting matter is subject to intense magnetic fields the QCD phase diagram boundaries 
are modified 0. Some of the most important changes concern the size and location of the first order chiral transition 
region since the results show that a strong magnetic field favors this type of transition. At the same time, at low 
temperatures, the value of the coexistence chemical potential decreases as B increases in accordance with the inverse 
magnetic catalysis phenomenon (ICM) Q. 

The low-T/high-/i region where a first order type transition is expected to occur is currently unavailable to lattice 
QCD evaluations (LQCD). However, the region high-T/low-/./, has already been exploited using LQCD simulations 
which indicate, in accordance with most model predictions, that the crossover observed at B = 0 persists when B / 0 
M- On the other hand, a major disagreement between recent LQCD results 0, [jj and model calculations regards 
the dependence of crossover pseudocritical temperature, T pc , on the strength B of the magnetic field. Specifically, 
the lattice results of Refs, p, [7[, performed with 2+1 quark flavors and physical pion mass values, predict an inverse 
catalysis, with T pc decreasing with B , while effective models predict an increase of Tpc with B. This problem has 
been recently addressed by different groups @,0 which basically agree that the different results stem from the fact 
that most effective models miss back reaction effects (the indirect interaction of gluons and B) as well as the QCD 
asymptotic freedom phenomenon. At the same time, other important aspects of the effects of strong magnetic fields 
on the QCD phase diagram have already been studied including the behavior of the coexistence chemical potential 
and the location of the critical end point (CEP) 0, the dependence of the CEP on strangeness, isospin and charge 
asymmetry m and also the internal structure of the phase diagram [0 . 

Regarding physical observables, the understanding of magnetized quark matter is particularly important at low 
densities and high temperatures, which is the relevant regime for the present heavy ion collisions experiments jldj ] , as 
well as at low temperatures and high densities, which is the regime concerning magnetars [ 13 ] . 

As far as heavy ion-collisions are concerned the presence of a strong magnetic field most certainly plays a role despite 
the fact that, in principle, the field intensity should decrease very rapidly lasting for about 1-2 frn/c only 0. The 
possibility that this short time interval may [0 or may not 0 be affected by conductivity remains under dispute. 
The effects of strong magnetic fields and their relation with the impact parameters have also been discussed 0 while 
the particle yields dependence on a constant external magnetic field has been investigated in a naive approach [0 . 
Another aspect related to the presence of strong magnetic fields at the early stages of the collisions is the anisotropy 
of photon production in heavy ion collisions at the RHIC energies 0 ] . The new PHENIX data brings some doubts on 
the conventional picture of thermalization and subsequent hydrodynamics or infers the possibility that a new photon 
production mechanism is possible. These works tell us that there is much to be done if a complete understanding of 
the effects of hadronic matter subject to strong magnetic fields is expected. 

When we look at the recent literature on magnetars, the controversy is already present at the level of calculat¬ 
ing the energy-momentum tensor. While some of the first works advocated that the pressure, having in mind the 
thermodynamical pressure obtained from the thermodynamical potential that relates pressure to density, should be 
isotropic f+ll+l , based on an interpretation given in Ref. 0] , other works were based on the fact that the energy- 
momentum tensor gives different contributions for the parallel and perpendicular pressure 0 - 0 ] (see Ref. 0 ] for 
a discussion based on LQCD). If two different pressures are indeed present in the system, the usual way of using the 
equation of state as input to the Tolman-Oppenheimer-Volkoff equations (TOY) [29j |. which determine the structure 
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of a spherically symmetric body of isotropic material in static gravitational equilibrium, to obtain compact objects 
macroscopic properties, as radii and masses, has to be done with care. What is normally done is to observe at what 
value of the magnetic fields, the two pressures start to deviate and then use the TOV equations up to this strength, 
so that in principle, the EoS used as input is practically isotropic [111 [53, [3G| • Another important aspect is related 
to the contribution of the electromagnetic interaction to the pressure(s) and energy density, a term proportional do 
B 2 , where B is the magnetic field strength. Since no field larger then 10 16 G has been observed at the surface of a 
magnetar, but according to the virial theorem one could expect fields as strong as ~ 10 18 G in the interior, an ad hoc 
exp onential density-dependent magnetic field was proposed in ref. 31] and widely adopted in subsequent works 0- 
l23l . I 32 III0I|. This ansatz, however, violates Maxwell equations. Another similar prescription for an energy density 
dependent magnetic field was proposed in Ref. [4l[ . 

To avoid the use of the TOV equations, the authors of Ref. [Hj treated the anisotropic pressure as a perturbation 
in a way similar to the Hartle-Thorn method, generally used for slowly rotating neutron stars. In a more complete 
treatment, the authors consider the anisotropy in solving Einstein’s field equations in axisymmetric regime with a fully 
general relativistic formalism [4^ - fi3 |. In both cases, once the macroscopic properties are obtained, a small increase 
of the maximum mass is found in contrast with the other previous works. Let us comment that the results computed 
from the LORENE modified code in Ref. [44|, generate a magnetic field in the central regions of the star of the order 
of 10 18 G, which seems reasonable but at the surface, its value is of the order of 10 17 G, which is too large in contrast 
with observational results that point to a maximum value of 10 15 G. 

According to classical books on gravitation [47],[48|, when anisotropies are present, the concept of pressure is not so 


well defined. Based on the concepts discussed in these two books, in Ref. [41| a small-scale chaotic field is used and 
the stress tensor is modified, so that the resulting EoS is also isotropic. A curious outcome is that the increase in the 
maximum stellar mass is also very small, as found in Ref. idil [iij . Hence, it is clear that there is no unique way of 
computing magnetic field effects on compact stars. 

An estimation of the maximum magnetic field intensity supported by a star before magnetic field stresses give rise 
to the formation of a black hole may be obtained equating the magnetic field energy of an uniform field in a sphere 
with the star radius R to the gravitational binding energy. A maximum field of the order of 10 18 G is obtained in 
agreement with the maximum fields obtained in the framework of a relativistic magneto-hydrostatic formalism, of the 
order of ~ 5 x 10 18 G with a nucleonic EOS [43[, or ~ 3 x 10 18 G in Ref. with an hyperonic EOS. It was suggested 
that a disordered field with ( B 2 ) > ( B ) 2 could possibly give rise to larger fields in still stable stars. The previous 
estimations referred to stars that are bound by gravitation. For self-bound stars larger fields could in principle exist 
[46j . Taking these numbers as indicative we next consider fields B < 1.5 x 10 19 G. 

In the present work we investigate, within the su(3) version of the Nambu- Jona-Lasinio |49j model the quark 
matter polarization and magnetization, the thermodynamical pressure, and the parallel and perpendicular pressure 
contributions obtained from the energy-momentum tensor. We consider both /3-equilibrium matter and quark matter 
with equal chemical potentials for the three flavors, which we call symmetric quark matter or matter with isochemical 
potentials throughout the text. The first scenario applies to neutron stars while the second is relevant to heavy ion 
collisions investigations. Some of these quantities are inputs for numerical codes that calculate the structure of neutron 
stars subject to strong magnetic fields. 

The paper is organized as follows: In Sections I and II, the general formalism and the resulting equation of state, 
developed in Refs. 0, 0] are revisited. In Section III, the expressions for the magnetization and the anisotropic 
pressures are shown, with some of the details given in Appendix A. In Section IV the results are shown and discussed 
and in Section IV the final conclusions are drawn. 


I. GENERAL FORMALISM 


In order to consider both symmetric quark matter and stellar quark matter in /3 equilibrium with strong magnetic 
fields we introduce the following Lagrangian density 

C = Cf+£ l -±F liU F'“', (1.1) 

which contains a quark sector, £/, a leptonic sector, Ci, and the electromagnetic contribution. The quark sector is 
described by the su(3) version of the Nambu-Jona-Lasinio model 

C f =-ijj f ['y IJ ,{id tl -qfA^-rhclijjf + C sym + C det . (1.2) 
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The C syrn and Cdet terms are given by: 

8 

C S ym = , ( 1 - 3 ) 

a =0 


G de t = -K {det/ [ip f (l + 75 )ip f ] + det f [ip f {l - 75 )V'/] } , (1-4) 

where ipf = (u,d,s) T represents a quark field with three flavors, rh c = diagj(?n K , rrid,m s ) with m u = 7 ^ m s is 

the corresponding (current) mass matrix while qj = diag(< 7 u , qd, q s ) is the matrix that represents the quark electric 
charges. In the same equation, Ao = -\/2/3 1 where / is the unit matrix in the three flavor space, and 0 < A a < 8 
denote the Gell-Mann matrices. The t’Hooft interaction term (Cdet) represents a determinant in flavor space which, 
for three flavor, gives a six-point interaction [50] and C syrn , which is symmetric under global XJ(Nf)i, x U(Nf)n 
transformations and corresponds to a four-point interaction in flavor space. The model is non renormalizable, and 
as a regularization scheme for the divergent ultraviolet integrals we use a sharp cut-off A in three-momentum space. 
The parameters of the model, A, the coupling constants G and K and the current quark masses, m u and m s , are 
determined by fitting /„■, m n , itik and m to their empirical values. We adopt the parametrization of the model 
proposed in [Hj: A = 631.4 MeV , m u = m d = 5.5 MeV, m s = 135.7 MeV, GA 2 = 1.835 and KA 5 = 9.29. 

The leptonic sector is described by 


Ci = ipi [q M (id 11 - qiA^) - mi] ^ , (1.5) 

where l = e,qi. One recognizes this sector as being represented by the usual QED type of Lagrangian density. As 
usual, Afj, and F yl , = dy,A v — d u A y are used to account for the external magnetic field. We are interested in a static 
and constant magnetic field in the z direction and hence, we choose the gauge A y = 


II. THE EOS 

We now need to evaluate the thermodynamical potential for the three flavor quark sector, Qf, which as usual can 
be written as 


n f = -Pf = £f-TS f -Y<L L fPf’ (2- 1 ) 

/ 

where Pf represents the pressure, £f the energy density, T the temperature, Sf the entropy density, [if the chemical 
potential, and pj the quark number density. A similar expression can be written for the leptonic sector. 

The total pressure for three flavor in f3 equilibrium is then given by 

= P f \ Mf + Pi\ mi ± — , ( 2 . 2 ) 

where our notation means that Pf is evaluated in terms of the quark effective mass, Mf, which is determined in a 
(nonperturbative) self consistent way while Pi is evaluated at the leptonic bare mass, m;. The term B 2 / 2, arises due 
to the kinetic term of the electromagnetic field, F llu F fJ ' v /4, in the original Lagrangian density. Within the formalism 
used in the present work, the sign of this term comes from the stress tensor and is shown in the section where we 
discuss anisotropy. Remark also that in the sequel our results will be presented in terms of (vacuum) subtracted 
pressures (AP) such that Pf = 0 at pf = 0 (/ = u,s,d) and Pi = 0 at /q = 0 (/ = e, p). With this normalization 
choice only the magnetic pressure, (±B 2 / 2), survives at vanishing chemical potentials. Again, a similar expression 
can be written for the leptonic sector, apart from the color index. The lepton masses are m e = 0.511 MeV and 
= 105.66 MeV. 

In the mean field approximation the pressure can be written as 

Pf = 9 U + Pd + — 2G(<(> 2 + <fid P (fig) + dK(j> u (j)d(t>s , (2-3) 

where the free gas type of term is 



d A p 


In ( -p 2 + M 2 f ) 


(2.4) 
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while the scalar condensates, (j)f are given by 


cj)f = (i/jft/jf) = -i 


d A p 

(27r) 4 (ji—Mf + ie) 


tr 


(2.5) 


According to standard Feynman rules for this model, all the traces are to be taken over color (N c = 3) and Dirac 
space, but not flavor. 

The effective quark masses can be obtained self consistently from 


Mi =rrii~- AG(f)i + 2 Kcj)j(j) k , 


( 2 . 6 ) 


with (i,j, k) being any permutation of ( u,d,s ). So, to determine the EOS for the su(3) NJL at finite density and in 
the presence of a magnetic field we need to know the condensates, </>/, as well as the contribution from the gas of 
quasiparticles, 9f. Both quantities, which are related by cf>f ~ ddf/dMf , have been evaluated with great detail in 
[2Cl l2l| so that here we just quote the results: 


0 f = ( 9) ac + 0™ a9 + 9f ed ^j 


M f 


where the vacuum contribution reads 


N P . 




(A + e A ) 
M f 


- e A A (A 2 + 4) | 


(2.7) 


( 2 . 8 ) 


and where we have also defined e A = A 2 + M 2 with A representing a non covariant ultraviolet cut off, the finite 
magnetic contribution is given by 


Q mag _ N c (\q f \B) 2 
f ~ 2tt 2 


1 X 2 

C'(—1 ,Xf) - ~(x 2 f - x f )hix f + '-j- 


(2.9) 


with Xf = M 2 / (2\qj\B) while (—1, a;/) = dC,{z,x f) / dz\ z —-\ where £(z,Xf) is the Riemann-Hurwitz zeta function. 
To take further derivatives, as well as for numerical purposes, it is useful to use the following representation for this 
quantity 


C ; (— x f) = C'(-IjO) + y[z/ - 1 - ln(27r) + ^ ( 2) (x f )\ , 


( 2 . 10 ) 


where ip( m \xf) is the m-th polygamma function and the Xf independent constant is £'(—1,0) = —1/12. The medium 
contribution can be written as 


Qmed _ T a j qf ^ c [ dp[\n(l+exp[-{E* f -p f )/T])+ln(l+exp[-{E* f +p, f )/T])] , 


k =0 


( 2 . 11 ) 


with ao = 1, afc>o = 2. In the above equation we have defined the energy dispersion 

E} = \JP 2 + s f(k, B) 2 , S f(k,B) = yjM} + 2\q f \Bk. 
When considering just the zero temperature case, Ea. (l2.1ip becomes: 


f) me d _ 'fT \9f\BN c 
a k- 

k—0 


47T 2 


Vfy/fif ~ Sf{k,B) 2 - s f (k, B) 2 !n 


df + V df - Sf{k,B) 2 
s f (k,B) 


( 2 . 12 ) 


where Sf(k, B) = M 2 + 2\qf\Bk. At T = 0, the upper Landau level (or the nearest integer) is defined by 

P 2 ~ M 2 P 2 ,f 


^ f,max — 


2\q f \B 2\q f \B- 


(2.13) 
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The condensates (f)f entering the quark pressure at finite density and in the presence of an external magnetic field 
can be written as 


where 


4>f = tt v f ac + <\>Y 9 + <t>T d ) 


Mf 


(2.14) 


(j) v f ac = - 


M f N c 

2n 2 


Aea — Mf In 


A + £a 
M 


f 


(2.15) 


jrnag M f \q f \BN c 

- 


27T 2 


lnr(x/) - iln(27r )+x f - ^ (2x f - l)ln(a:/) 


(2.16) 


and 


,med V, M f \q f \BN c [ + °° (/+ + /-) 

<t>f = 2^ ak -139- / dp- 


k—0 


47T 2 


E f 


(2.17) 


where the Fermi distribution functions are 


The quark density reads 


f± = 1/(1 + exp[(^/ T Pf)/T}} . 


k v^‘ \qf\BN c [ +oc , (f ,, 

Pf= 2^ ak 2n 2 ' d p(f+-f-). 


k—0 


At T = 0, Eas. (12.17l) and (12.1911 become: 


(2.18) 


(2.19) 




imed \ ^ 

h = 2^ 


CCfc- 


M f \q f \BN c 


k-0 


2ir 2 


111 


Pf + \Jp f ~Sf(k,B) 2 
s f(k, B) 


( 2 . 20 ) 


and 


«= 'f 


k—0 


where fcpj = ^ P 2 f — s f{k, B) 2 . The entropy density Sf = —(d£l/dT) is 


*/ = -£ £ 


Oik 


} fc=0 


\q f \BN c 

Air 2 


dp lf + In (/+) + (1 - U) In (1 - /+) + (/+ O /_)] . 


( 2 . 21 ) 


( 2 . 22 ) 


The corresponding leptonic contributions can be trivially obtained from the above quantities by replacing Mf —> mi , 
19/1 I®I an d Pf ► Pi- Since the leptonic masses are unaffected by strong interactions one considers their bare 

values which do not depend on T, p and B as opposed to the effective Mf masses. Therefore, the only piece which 
effectively contributes to the subtracted pressure, defined as A P = P(T, p, B) — P( 0,0, B ), is 


k f,‘ 


pmed _ rj -r ^ ^ 


Oik' 


k—0 


\ qi\B 

Air 2 


r +oo 


dp [In (1 + exp [-{Ei - pi)/T ]) + In (1 + exp [-{Ei + pi)/T])\ , 


(2.23) 


from which the leptonic density can be written as 


^ f,max 


pi= a k 


\m\B 


/*+oo 


2t r 2 / 

k =0 d-oo 


dp{l+ -1-) , 


(2.24) 
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where m represents the leptonic chemical potential. The quantities Ei and l± can be obtained from their quark 
counterparts by the replacements already mentioned. 

At vanishing temperatures the above expressions become: 


r>med 


fJ- fcl,max 

= E E 

l=e k=0 




\<n\ B 

47T 2 


Pi 


n n m2, ( « + - si(k,By 

B) 1»(-- 


and 


Pi 


E ak 


k —0 


\qi\ B 

2n 2 


k F ,i(k,si ) , 


(2.25) 


where k F ,i{k,Si) = rj[if - Si(k,B) 2 . 


III. THE ANISOTROPY IN THE PRESSURE 


The parallel and the perpendicular components of the pressure can be written in terms of the magnetization, 
M = <9A P/dB, as 0,0,0: 


d2 d2 

P|l = AP - — and P ± = AP- MB + — (3.1) 

where AP stands for the already defined subtracted pressure. For a magnetic field in the z direction, the stress 
tensor has the form: diag(P 2 /2, P 2 /2, —P 2 /2) and this explains the difference in sign appearing in the parallel and 
perpendicular pressures. For the leptonic sector one easily gets 


which, at T = 0 becomes: 


d prned _ prned |^|£ 

dB ~ ~B 47T 2 

d pmed _ prned £| ?J | 

dB ~ ~B 2^ 


k f^max 

E a fc( fc i®i) 

k—0 



dp -j^[ l + + 1 -} 


k-max 

Y, OL k In 


k =0 




( fc M). 


(3.2) 


(3.3) 


whereas for the quark sector one obtains 
dPt 

= O'u + ®'d + 9's ~ 4 G((f) u 4>' u + <j>d<t>d + ^sfis) + 4 K(tj) u <\) d <\) a + <j) u (j)' d (j) s + (frufidfis) , (3.4) 

where 

e' f = (e' f vac + e' f ma9 + e' f med ) Mf , ( 3 . 5 ) 

and 

4>'f = {<t>' f vac + t'™ 9 + 0/ med K ■ (3.6) 

The primes denote derivatives with respect to B and the explicit form of each term can be found in the appendix A. 


IV. RESULTS AND DISCUSSION 

In this section we present and discuss results concerning several properties of /3-equilibrium quark matter and 
symmetric quark matter (with equal chemical potentials). We first discuss some general properties of quark matter, 
in particular, its particle content and particle polarization. We then investigate how the magnetization of quark matter 
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FIG. 1: /3-equilibrium quark matter (left panels) and symmetric quark matter (right panels): a) and b) particle 
polarizations for several densities defined in terms of the saturation density po = 0.16 fm , c) and d) particle 
densities, absolute in c) and relative in d), for 10 17 G, 10 18 G and 10 19 G magnetic field intensities, respectively, from 
the thinner to the thicker lines; e) total polarization onset density for the s quark and the electron in ^-equilibrium 
matter and f) for the u and d quarks in symmetric matter. For the densities representing u and d quarks the total 
polarization occurs for stronger fields and in symmetric matter the s quark sets in at larger densities. 



























changes with the magnetic field intensity and finally we discuss the parallel and perpendicular pressure contributions 
obtained from the energy-momentum tensor. 

In Fig. [l] it is shown how the magnetic field affects /3-equilibrium matter (left panels) and symmetric quark matter 
(right panels), in particular, the quark and electron polarizations, their densities and onset of total polarization. The 
^-equilibrium matter onset of the s quark occurs just below 0.7 fm” 3 , see Fig. [Tc] and, therefore, close to these 
densities the s quark density is small and feels strongly the magnetic field, total polarization being attained with 
B < 10 18 G. In denser matter, the larger densities of s quarks require larger magnetic field intensities for a total 
polarization. In symmetric quark matter the s quarks set in above 0.9 fm -3 , a density that presently is not attained 
in the laboratory. In /3-equilibrium matter u and d quarks are not totally polarized for fields below 10 19 G, however in 
symmetric quark matter u and d quark total polarization occurs at small densities that do not exist in quark stellar 
matter with a surface baryonic density which is of the order of 0.3 fm -3 . In Fig. [la] we also show information on 
the electron polarization. The magnetic field increases the electron content and for 10 18 G their density is practically 
constant and equal to ~ 0.01 fm” 3 , see Fig. [Tc] At the onset of the s quark the density of electron has always 
a maximum, above which the electron fraction decreases. This means that electrons are totally polarized for fields 
B > 9 x 10 17 G. 

For totally polarized matter, all particles lie on the lowest Landau level (LLL). In this case, the dependence of 
the pressure on the magnetic field intensity of a gas of free particles occurs only through a multiplicative factor that 
defines the LLL degeneracy, and the magnetization is independent of B. For the NJL model the interaction terms 
give a nonlinear dependence to the pressure even above for totally polarized matter. This contribution becomes small 
when the chiral symmetry is restored. 

In order to determine the two pressure contributions, the magnetization is obtained using Eq. (16.11) . In Fig. [2] 
this quantity is plotted for /3-equilibrium matter, Figs. I2al and Eel and quark matter with equal chemical potentials, 
Figs. I2bl and [2dl For each case we have considered a set of densities of interest: a) the surface baryonic density of a 
quark star is ~ 2po and in the interior we may have densities larger than 5po; b) in heavy ion collisions we may have 
densities below po and do not expect densities as large as 5po- The magnetization has a term that explodes whenever 
Pf = Sf, i.e., whenever B approaches an/0 Landau level from below, see Eq. (16.71) and the corresponding discussion 
in Ref. [25j. The spikes in this figure occur precisely at this values of B. A small number of spikes occurs if only a 
few LL are occupied, e.g. if the magnetic field is very strong with respect to the Fermi momentum of the particle. In 
Figs. Ha] and (20 for the densities represented this occurs for the larger fields. For the symmetric matter at smaller 
densities, fewer spikes are obtained for larger magnetic fields. It is interesting to notice that for /3-equilibrium matter 
there is a range of densities between 0.5 and 0.7 fm -3 where for all the fields shown the number of oscillations is 
small: this occurs before the onset of the s quark and after the u and q quark restoration of chiral symmetry. In Fig. 
Eel for densities lower than 0.3 fm” 3 and in between 0.8 and 1.1 fm” 3 , one can clearly see that oscillations with high 
frequencies are modulated by smaller frequency oscillations. The same is seen in Fig. I2dl for densities larger than 
0.3 fm” 3 . This superposition of fluctuations with different frequencies is due to a mixing of particles with different 
charges and masses. In /3-equilibrium matter the oscillations defined by the u quark are wider apart because the u 
quark density is smaller and its charge larger, both effects adding to a reduction of the number of occupied LL, see 
Fig. HB This same effect if also present in Fig. [2d] where, for similar u and d quark densities, the charge of the u 
explains the difference. 

In Fig. [3] the two pressure contributions defined Eq. (13.11) . except for the terms proportional to B 2 , are plotted for 
both /3-equilibrium matter (left panels) and symmetric matter (right panels). As in the previous discussion for each 
scenario we take the same representative values of the magnetic field intensity or baryonic density. 

For fields B < 10 18 G both contributions, the parallel and the perpendicular pressures, are practically coincident. 
However, due to the magnetization contribution entering the perpendicular pressure, discontinuities occur whenever 
a new LL level is reached, giving rise to a series of discontinuities that would correspond to unstable regions, since 
dJAfdB < 0. Let us point out that in the region 0.5 < p < 0.7 fm” 3 , and for B < 10 18 G, these discontinuous con¬ 
tributions to the perpendicular pressure contribution are practically zero, because, as seen before, the magnetization 
is very small in this range of densities. As compared with calculations performed with the MIT bag model [30j ]. the 
behavior is the same, except that the deviation of the two pressures takes place at even larger magnetic fields with the 
MIT bag model. Nevertheless, had we considered not only the matter contribution but also the contributions from 
the the terms proportional to B 2 , the deviation would start at lower magnetic field intensities. 

One should bear in mind that the discontinuities seen in Fig. [3] are washed out at finite temperatures as the ones 
that we expect in heavy ion collisions [52|. As suggested in Refs. [U, [25], |46| the unstable regions may give rise 
to domains with disordered fields that would allow to continuous phase transitions between regions with a different 
number of LL. 
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(c) (d) 

FIG. 2: Magnetization for a) and c) /3-equilibrium matter and b) and d) isochemical potential quark matter. 
Densities of interest for each scenario are considered. In a) and b) the magnetization is shown as a function of B 

and in c) and d) as a function of the baryonic density. 


V. FINAL REMARKS 

In this work we have examined the effects of strong magnetic fields in quark matter, as described by the NJL model, 
paying special to effects due to pressure anisotropy. Two scenarios were investigated: matter in /3-equilibrium, as the 
one possibly present in quark stars or in the core of hybrid stars and symmetric matter with isochemical potentials, 
as the one present in heavy ion collisions. For the first case, large densities and magnetic fields up to the order of 
10 18 G are of interest whereas for symmetric matter, the relevant densities are somewhat smaller while the magnetic 
fields can be larger. 

Part of our work has been devoted to the analysis of quark matter constituents, the onset of s quarks and how much 
polarized matter can be present when different field intensities are considered since these are important ingredients 
for stellar modeling. 

Magnetization, which is an important quantity for the description of magnetized matter, has also been considered 
in the present work and by investigating this quantity we conclude that the amount of spikes, related to the filling of 
the LL, depends quite substantially on the density and on the scenario examined. Of particular interest is the fact 
that very few spikes are seen for densities between 0.5 and 0.7 fm~ 3 in /3-equilibrium matter, no matter the intensity 
of the magnetic field, because the particles that constitute matter occupy only a few LL. For symmetric matter, the 
same pattern is found for densities larger than 0.5 fm , when the number of occupied LL are small for the three 
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(a) 



P equil 



(b) 



Log 10 [B(G)] 


(c) (d) 

FIG. 3: Parallel (lines) and perpendicular (marks) pressures for a) and c) /3-equilibrium matter and b) and d) 
isochemical potential quark matter. Densities of interest for each scenario are considered. In a) and b) the pressures 
are shown as a function of the baryonic density and in c) and d) as a function of B . For fields B < 10 18 G the 
perpendicular press falls on top of the parallel contribution except for the van Alphen oscillations due to the 

magnetization contribution. 


different quarks. 

This effect in /3-equilibrium matter is due to the late onset of que s quark, precisely for p > 0.7fm~ 3 . Due to its 
large constituent mass magnetic fields satisfying B < 10 18 G, as the ones possibly existing inside magnetars, give rise 
to the filling of many LL and large contributions to the magnetization and perpendicular pressure, in contrast to the 
behavior at densities just below the s quark onset. This effect is clearly seen comparing the pressure for p = 3po 
corresponding to u and d quark matter with restored chiral symmetry with the pressure obtained for either p = 2po 
or p = 5po as a function of B: in the first case chiral symmetry is still not completely restored and in the second the 
onset of s quark has occured. If a quark model with complete chiral symmetry restored such as the MIT bag model, 
or a model without the strangeness degree of freedom such as the su(2) NJL is used to describe /3-equilibrium quark 
matter, a different behavior will probabily occur for B < 10 18 G and p > 0.5 fm -3 right until the star center, e. g. 
p ~ 1.2 fm -3 [2ll [53| , in particular, both pressure contributions will be coincident. 

Therefore within the su(3) NJL model, a non negligible effect of the magnetic field on the quark star structure close 
to the surface and in the interior is expected even for B ~ 10 18 G. Taking as reference the calculation done in [54j 
using pure toroidal magnetic field equilibrium models of relativistic stars for both non-rotating and rotating hadronic 
stars, fields as large as 10 18 G were obtained on the equatorial plane deep inside the star. In this calculation, however, 
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the magnetic field effects on the EOS have been neglected. 

We have then looked at the pressure anisotropy obtained from the components of the energy-momentum tensor 
that define the parallel and the perpendicular pressure contributions and examined the relation between parallel and 
perpendicular pressures for both scenarios described above when the magnetic field is fixed and also when the baryonic 
density is kept constant. We have observed that the larger the magnetic field intensity, the larger the discontinuities in 
the perpendicular pressure. It is important to note that the magnetization is not responsible for the complete picture, 
because it comes multiplied by the magnetic field in the calculation of the perpendicular pressure. In agreement 
with Refs. [26j, j30| we have observed that when the densities are fixed, the parallel and perpendicular pressures are 
practically coincident up to very large magnetic fields, what could justify the use of isotropic matter hydrostatic 
equations in stellar calculations. However, the integration of the full relativistic hydrostatic equations still requires 
the inclusion of the magnetic field contributions which involve a B 2 term that gives rise to a quite large effect. 
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VI. APPENDIX A - DERIVATIVES WITH RESPECT TO B 

As already shown in the text, the magnetization is given by: 


but in equilibrium 


then 



dn \ 

dW s ) 




= 0 



( 6 . 1 ) 


Thus, in the following expressions we do not need to include terms containing which also means that the vacuum 
contributions d^ vac and (j>'^ vac trivially vanish. Then, 
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For T = 0 the above relation becomes: 
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A straightforward evaluation yields 
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(6.5) 


where ip°(xf) = is the digamma function. The in medium contribution reads 
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which, at T = 0, can be written as 
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